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We study the scaling behavior of the step scaling function for SU(3) gauge theory, employing 
the renormalization-group improved Iwasaki gauge action and the perturbatively improved Liischer- 
Weisz gauge action. We confirm that the step scaling functions from the improved gauge actions 
■ agree with that previously obtained from the plaquette action within errors in the continuum limit 

' at both weak and strong coupling regions. We also investigate how different choices of boundary 

counter terms for the improved gauge actions affect the scaling behavior. In the extrapolation to the 
continuum limit, we observe that the cut off dependence becomes moderate for the Iwasaki action, 
^ \ if a perturbative reduction of scaling violations is applied to the simulation results. We also measure 

the low energy scale ratio with the Iwasaki action, and confirm its universality. 



I. INTRODUCTION 
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The strong coupling constant is one of the fundamental parameters of QCD. The current world average leads to 
a Ms( TO - z ) = 0.1172(20) [1]. Lattice QCD calculations have a potential ability to determine the strong coupling constant 
from an experimental input at low energy scales. In practice, however, one must relate the high energy perturbative 
QCD scale to the low energy hadronic scale. Alpha Collaboration proposed the Schrodinger functional (SF) scheme 
as a vehicle for this purpose [2,3], and it has been successfully applied to lattice QCD in various aspects [4-7]. One 
of the most recent results related to our study is the running coupling constant of two massless flavor QCD reported 
• • ! in Ref. [8,9] 

Recently CP-PACS and JLQCD Collaborations have started a project for Nf = 3 QCD simulations [10-14]. These 
simulations are essential to understand the low energy QCD dynamics for the real world in which three light quarks 
exist. One of the targets of the project is to evaluate the strong coupling constant a-^r in Nf = 3 QCD using the 
SF scheme. In the project Iwasaki gauge action [15] is employed to avoid the strong lattice artifacts of the plaquette 
gauge action found in Nf = 3 simulations [11]. 

In a previous study [16], as our first step toward evaluation of ot^r for Nf = 3, O(a) boundary improvement 
coefficients in the SF scheme have been determined for various improved gauge actions up to one-loop order in 
perturbation theory. In addition the scaling violation in the step scaling function(SSF) for the coupling have been 
analyzed perturbatively. In the present paper, as the next step, we investigate the lattice cut off dependence of the 
SSF non-perturbatively in quenched lattice QCD simulations with improved gauge actions. The renormalization- 
group improved Iwasaki gauge action and the perturbatively improved Luscher-Weisz gauge action are employed. We 
investigate the effect of various choices for boundary improvement coefficients in detail, to find the best choice, which 
will be used in our unquenched simulations in the future. We also confirm the universality of the SSF and the low 
energy scale ratio, by comparing our results with the previous ones obtained by ALPHA Collaboration [17,18] 

The rest of this paper is organized as follows. In Sec. II, after a brief introduction of the SF scheme and its 
extension to improved gauge actions, we specify the action and the O(a) boundary improvement coefficients used in 
our simulations. We then define the Schrodinger functional coupling constant, the step scaling function and the low 
energy scale ratio. In Sec. Ill, we give details of simulations and present our results with improved gauge actions 
for various choices for O(a) improvement. In Sec. IV, we investigate the lattice cut off dependence of the step 
scaling function and the low energy scale ratio, and carefully take the continuum limit of these quantities, in order to 
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confirm their universality. Our conclusion is given in the last section, together with a discussion toward Nf = 2 and 
3 simulations. 



II. PRELIMINARIES 
A. Schrodinger functional 

The SF scheme introduced by ALPHA Collaboration is a powerful tool to probe the energy evolution of physical 
quantities. In the SF scheme, the theory is defined on a finite box of size L 3 xT with the periodic boundary condition 
in the spatial directions and the Dirichlet boundary condition in the time direction. We set T = L throughout this 
paper. In the pure SU(3) gauge theory with Wilson plaquette action S[U], the Schrodinger functional is given by 




(11.1) 



where the link variables U(fi,x) for the gauge fields satisfy the boundary conditions 

U(x,k)\ Xo=0 = exp{aC}, U(x 1 k)\ Xo=L =cxp{aC'}. (11.2) 

Here a is lattice spacing, and C, C are spatially constant diagonal matrices, which depend on the background field 
parameters r\ and v [19]. 

An extension of the SF scheme to the improved gauge actions was first discussed by Klassen [20] in terms of 
a transfer matrix construction [21]. In this formulation, each boundary consists of two time slices, to achieve the 
tree-level 0(a 2 ) improvement. 

In this paper, however, we adopt an alternative formulation [22], which achieves the tree-level 0(a) improvement 
with only one time slice at each boundary. The dynamical variables to be integrated over are independent of the form 
of the action, whether plaquette or improved, and consist of the spatial link variables U (k, x) with xq — a, ■ ■ ■ , L — a 
and temporal link variables U(0, x) with x a = 0, ■ • • , L — a on the cylinder with volume L 3 x L. This formulation is 
implemented more easily in numerical simulations. 



B. Gauge action and O(a) boundary improvement coefficients 

The improved action we employ includes the plaquette and rectangle loops, and is given by 

SitnAU] 4E W (C,g 2 )2C(C) + \ ]T Wi(C, g 2 )2£(C), (II.3) 

with 

C{C) = ReTr[/-C/(C)], (II.4) 

where Wi is a weight factor to be specified later and U(C) is an ordered product of the link variables along a loop 
C contained in a set So (plaquette) or Si (rectangular) . Sq and Si consist of all loops of the given shape which can 
be drawn on the cylindrical lattice with the volume L 3 x L. The loops involve the "dynamical links" in the sense 
specified above, and spatial links on the boundaries at x — and x — L. In particular, rectangles protruding from 
the boundary of the cylinder are not included. 

One needs to choose the weight factors appropriately to achieve the one- loop level 0(a) improvement. Among 
various possible choices, ours is given as follows. 



W (C,g 2 ) 



c c s ((/q) for C G P s : Set of plaquettes that lie completely 

on one of the boundaries, 
Cocf (<7q) for C G Pt : Set of plaquettes that just touch one (H.5) 

of the boundaries, 
cq for C € Pother : otherwise, 
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with 



) for C G R s : Set of rectangles that lie completely 

on one of the boundaries, 
W\{C,gl) — { Ciqf(<7o) f° r C e i?^ : Set of rectangles that have exactly two (II. 6) 

links on a boundary, 
ci for C € -Rother : otherwise, 

c Ct P (5o 2 ) = Cb(l + cf (1) 5o 2 + O(. 9o 4 )), (II.7) 
Cl cf ( 5o 2 ) = Cl (3/2 + cf (1) . 9o 2 + O( 3o 4 )), (II.8) 

where the coefficients c and ci of the improved gauge action are normalized such that Co + 8ci = 1. In this paper 
we consider not only the Iwasaki action (ci = —0.331, C2 = C3 = 0) [15] but also the Liischer-Weisz (LW) action 
(ci = — 1/12, C2 = C3 = 0) [23] for comparison, since our perturbative analysis [16] shows that the LW action has 
a fairly small lattice artifact in the step scaling function. We call cf (g 2 ) and cf(g 2 ) 0(a) boundary improvement 
coefficients. The assignments at the t = boundary are shown in Fig. 1. The leading term of the 0(a) boundary 
improvement coefficients in eqs.(II.7) and(II.8) can uniquely be determined from two requirements that the tree-level 
O(a) improvement is achieved and the lattice background field satisfies the equation of motion at the boundaries [22]. 
On the other hand, for the one-loop boundary terms, we can freely set a relation between c t v ; and c t , since there 
is only one requirement for the one- loop 0(a) improvement. 

Let us see how we specify the one-loop boundary terms. In Ref. [16], one finds the following relation to achieve the 
one-loop 0(a) improvement: 

c cf (1) + 4c lC f (1) =A 1 /2, (11.9) 

where A\ is the coefficient of the a/L term in the one- loop correction (L/a) to the SF coupling. In our simulations 
we consider two choices: one called condition A is given by 

cf«= 2c r«, (ii.io) 

and the other called condition B is specified by 

cf (1) =0. (11.11) 

The difference between the conditions A and B is an 0(a 5 ) contribution in the one- loop correction to the SF coupling 
[16]. Although this difference is tiny at one-loop order, it may become larger at higher orders. The values of one-loop 
boundary terms for each condition and A\ are given in Table I. For the LW action, the difference between the two 
conditions is small, so wc do not carry out simulations with the condition B. 

C. Schrodinger functional coupling 

The SF with the improved gauge action is given by 

Z = er r = J D[U]e- Sl ^ [u] , (11.12) 

where we impose the same boundary condition eq. (II. 2) for the link variables as in the case of the Wilson plaquctte 
action. The SF coupling is defined through the free energy T in eq. (11.12) 



'AC 

- g l F (L) = k/r>\ v=v=Q = k/{ — 



(11.13) 

•q— v— 



L N 2 



where k is a normalization constant 

fc=12^J [co (sin27 + sin7) + 4ci (sin 47 + sin 27)] , (II. 14) 

Hz)VS)- ou« 

The renormalized coupling <?g F (L) depends only on the scale determined by the box size L. The //-derivative of the 
loop touching the boundaries (the left-most three loops in Fig. 1, for example), contributes to the observable dS/drj 
ineq. (11.13). 
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D. Step scaling function and Low energy scale ratio 



The step scaling function (SSF) describes the evolution of the renormalized coupling under a finite rescaling factor 
s (we take s — 2 in the following) 

a(2,u) = f{2L)\_ 2(L)=u . (11.16) 

By choosing the n + 1-th initial value of a(2, u n+ \) such that u n+ \ = cr(2, u n ), the non-perturbative evolution of the 
running coupling can be constructed successively in order to cover a wide range of the energy scale. 

The SSF <t(2, u) in the continuum theory is obtained by the continuum limit of a lattice SSF S(2, u, a/L) 

a(2,u) = lim E(2,u,a/L). (11.17) 

a/L^O 

In this paper, we study the SSF at a weak coupling u — 0.9944 and a strong coupling u = 2.4484, where our results 
can be compared with those of ALPHA Collaboration. 

To fix the scale in a physical unit, one needs to relate the box size L prescribed at a certain value §s F (L) to some 
reference scale. Following the conventional way, we set L — L max defined implicitly 

5i F (£max) = 3.480, (11.18) 

and adopt Sommer's scale tq [24] as the reference scale. Eventually this amounts to computing the ratio L max /ro and 
extrapolating it to the continuum limit 

imax/ro = lim (L max /a) x (a/r ). (11.19) 

0/Lmax^0 



III. SIMULATION DETAILS AND RESULTS 



We follow the calculation procedure of Ref. [19]. Simulations for the SSF on larger lattices are performed on 
CP-PACS using 4 partitions of 64 PU's, while r /a are calculated with 2 partitions of 512 PU's. 

As mentioned in Sec. II C, the SF coupling is obtained by calculating the observable dS/drj for gauge configurations 
with the SF boundary condition. When the number of spatial lattice points, L/a, is a multiple of 4, the gauge 
configurations are generated by a combined five-hit pscudo-heat-bath algorithm and an over-relaxation algorithm 
(HB). The combination of one pscudo-heat-bath update sweep followed by Nor = L/2a over- relaxation sweeps is 
called an iteration. The measurement is implemented after each sweep, i.e. (1 + L/2a) measurements are made per 
one iteration. Because of a restriction of the HB method optimized for the improved gauge action on CP-PACS, we 
employ the hybrid Monte Carlo (HMC) algorithm for L/a being different from multiples of 4, for instance L/a = 6. 
The step size for the molecular dynamics is adjusted to achieve an acceptance rate in a range from 0.7 to 0.8. The 
measurement is made for every trajectory. 

Our computations for the renormalized coupling g~s F (L) are carried out on lattices L/a — 4, 6, 8, 12. In this 
calculation, a re- weighting technique is used for a tuning of [25] such that <7g F (L) becomes a certain prescribed 
value u for each L/a. And then, using the same 0, a computation on a lattice with twice the linear size 2L/a gives 
g<j F (2L). The results are summarized in Table II and III for the weak(u = 0.9944) and strong(w = 2.4484) couplings, 
respectively. Errors in both g$ F (L) and cjg F (2L) are estimated by a jack-knife method. The bin size for jack-knife 
errors is 100 iterations for the HB and 500 trajectories for HMC, respectively. The precision in g^ F (L) is attained 
by accumulating around 120000 — 140000 iterations on the lattices L/a — 4, 8, 12, and 300000 trajectories on the 
lattices L/a — 6. As for g~g F (2L), the number of iterations is around 40000 — 80000 to achieve their precision. Errors 
in g<j F (L) are propagated into S(2, u, a/L), the lattice SSF, where u is the central value of <?g F (L). A formulae of the 
error propagation using a perturbative expansion of the SSF can be found in Ref. [26] . 

We performed an additional set of simulations with the Iwasaki action to determine the low energy scale. The tuning 
of to the conventional point gg F (L max ) = 3.480 and the error analysis are made in the same way as mentioned above. 
In Table IV we list the results, which will be used in Sec. IV B as the first factor on the right hand side of eq. (11.19). 
To complete the scale determination, one needs the second factor in eq. (11.19). In addition to the previous results of 
r /a [27-29], we carried out simulations at = 3.00 and 3.53 to cover the range of in Table IV. Analysis procedures 
for the static quark potential and extraction of r /a parallel those in Ref. [29] . The simulation parameters and results 
in this work are shown in Table V. To avoid finite size effects, we followed a criterion [18] that the parameter and 
L/a are chosen such that L/ro ~ 3.3. Following the above reference, the number of over-relaxation sweeps are taken 
to satisfy Nqr ~ l-5(r /a). 
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IV. CONTINUUM EXTRAPOLATION 



A. Step scaling function 



In this subsection we investigate the cut off dependence of the SSF and perform the continuum extrapolation 

a(2,u)= lim Y,(2,u,a/L). (IV. 1) 

a/L^O 

The lattice SSF S(2,«, a/L) as a function of a/L at the weak coupling u — 0.9944 is shown in Fig. 2 (a) and (c) 
for the Iwasaki action and LW action, respectively. For the Iwasaki action, even after the one- loop 0(a) improvement 
with either the condition A or B, the scaling violation is still rather large, which makes the extrapolation to the 
continuum limit difficult. To improve the scaling behaviour of the SSF, we apply a perturbative removal of the lattice 
artifacts suggested in Ref. [30] given by 

where Y.^(2,u,a/L) is the SSF (simulation raw data) with the "fc" -level 0(a) improvement coefficient (e.g., k = 0: 

tree level 0(a) improvement case, k = IA: one-loop 0(a) improvement with condition A case, etc.). S[ k \a/L) is the 
one-loop relative deviation, given by 

^)(2, U ,a/L) = S(fc) ^ U ^-^ U ) 

= S {k) (2,a/L)u + 0(au 2 ), (IV.3) 
whose numerical values are given in Table VII. This method eliminates not only O(a) but also 0(a n ) with n > 1 

(k) 

lattice artifacts at one- loop order. Figure 2 (b) shows the cut off dependence of T,\ (2,u, a/L). Indeed the scaling 
violations are much reduced by this method, so that we can reliably take the continuum extrapolation linearly in a as 

Z[ k) (2, u, a/L) = ct« (2, u) + u[ k) (u)a/L, (IV.4) 

where a^(2, u) and uj[ k ^ (u) are fit parameters. In Table VIII we quote the extrapolated value for the Iwasaki action, 
which is obtained by a simultaneous fit for k = 0, 1A and IB data with the constraint that they agree in the continuum 
limit a i - k) (2,u) = a(2,u). The fit has a good x 2 / N dof as listed in Table VIII. 

As shown in Fig. 2 (c), the scaling violations arc quite small for the LW action. This is consistent with the fact 
found in Ref. [16] that the lattice artifacts arc quite small at one loop. Moreover, the difference between the tree 
level and one- loop 0(a) improvement is invisible in this precision, as a result of the smallness of the improvement 
coefficients cf ^ and As shown in Fig. 2 (d), the perturbative removal of lattice artifacts has almost no effect 

except for L/a = 4, since 6[ k \2,a/L) with L/a = 6, 8, 12 is quite small. In Table VIII we quote the extrapolated 
value for the LW action, obtained by a linear fit to data of the one-loop 0(a) improved action with the perturbative 
removal of lattice artifacts. For comparison results of ALPHA Collaboration are also included [17]. 

Results at the strong coupling u — 2.4484 are plotted in Fig. 3. For the Iwasaki action, the one-loop 0(a) 
improvement shows large lattice artifacts for both conditions A and B, particularly for the coarse lattice (see Fig. 3 
(a)). As shown in Fig. 3 (b), the perturbative removal of lattice artifacts well reduces the scaling violation in the case 
of the condition B, but it still remains rather large for the condition A. Therefore, we include a quadratic term in the 
fitting form, 

X {k ) (2, u, a/L) = a (fe) (2, u) + J k) (u)a/L + J k) {u){a/L) 2 , (IV.5) 

for the data of k = 1A and IB, while we use the linear fitting form eq. (IV.4) for the data of k = 0. The extrapolated 
values obtained with the constraint a^(2,u) = a(2,u) for a unique continuum value is listed in Table VIII. We 
note that \uj < 2 A \u)/uj^ A \u)\ s=s 0(10) in the fit for the condition A. This suggests that the condition A accidentally 
enhances the coefficient of 0(a 2 ) term. It does not necessarily mean, however, that the one-loop 0(a) improvement 
itself is inefficient at this coupling constant. Indeed the one- loop 0(a) improvement with the condition B shows good 
scaling behavior. As for the LW action, Fig. 3 (c) and (d) show that neither the one- loop 0(a) improvement or the 
perturbative removal works effectively. Concerning the extrapolation we simply use the same procedure as in the 
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weak coupling case, i.e. the linear fitting form to the perturbative removal data for one-loop 0(a) improvement. The 
result is given in Table VIII. 

We observe in Table VIII that the three values obtained with the Iwasaki and LW action in the present work and 
that of ALPHA Collaboration [17] at the weak coupling are consistent within la. At the strong coupling, the value 
for the LW action undershoots relative to the others by 1.5-2ct. We think that the latter disagreement is caused by a 
large lattice artifact for the LW action, which makes the choice of the fitting form difficult. For example, if we assume 
that O(a) errors for the LW action are negligible, we can obtain a result consistent with the values of the other actions 
within la, by using a purely quadratic fitting form. Further investigation are needed to clarify this point. 

B. Low energy scale ratio 

We now combine L max /a and a/ro for the Iwasaki action to form the ratio I ma x/fo, and extrapolate it to the 
continuum limit: 

L max /r Q = um (L max /a) x (a/r ). (IV.6) 

«/im>x^O 

In the fourth column of Table IX, we give the first factor, which is taken from Table IV; the error of L max /a is 
estimated by propagating that of ga,-p(L max ). 

The second factor r /a is given in the third column of Table IX. This is obtained by an interpolation of the results 
for r /a in Table VI using a polynomial [31] 

In(o/r ) = ci + c 2 (/3 - 3) + c 3 (/3 - 3) 2 . (IV.7) 

The fit, plotted in Fig. 4, gives 

ci = -2.193(6), c 2 = -1.344(7), c 3 = 0.191(24), (IV.8) 

with x 2 /A r dof = 4.10/6 in the range 2.456 < fi < 3.53. The error of r /a in Table IX includes both statistical and 
systematic ones. We take the central value of r /a from the result of the fit above, and estimate the systematic error 
from the difference of r$ /a between the central value and the result of another fit including a C4 (f3 - 3) 3 term. 

The combination of the two factors for various (3 values are listed in the fifth column in Table IX. The /3 dependence 
of L max /r can be considered as lattice cut off effects. For an extrapolation to the continuum limit, we use a fit form 

a/r = ha/L max + b 2 {a/L max f + o 3 (a/L max ) 3 , (IV.9) 

where bi(i — 1,2,3) are fit parameters and 61 is the continuum value of the low energy scale ratio, rather than fitting 
L max /r as a function of a/r . In this way one can avoid the correlation of errors which complicates the latter fit 1 . 
We apply eq.(IV.9) to three sets of data, i.e., data for the tree level 0(a) improvement and those of the one- loop 0(a) 
improvement with the conditions A and B. For the first set of data, we set o 3 = (a linear fit) and exclude the point 
at a/L = 1/4. An alternative fit including that point and allowing a non-zero 03 yields a consistent value for b\ within 
errors. However, we think that the fit is not so reliable since xV-^dof is too small and the linear terms are rather 
large. Therefore we exclude the point and use the linear fit for the remaining 3 points. In Fig. 5 L max /r is plotted 
as a function of a/L max . Dashed lines are the fit curves eq.(IV.9) divided by a/L max , and the point at a/L max = 
shows the extrapolated value b\. 

The extrapolated value is given in Table X, together with the previous result for the standard Wilson plaquette 
action [18]. While rather large lattice artifacts are observed for both standard Wilson plaquette action and Iwasaki 
action, the extrapolated values agree within errors. 



1 Since the error on L max /a is small, one may perform a continuum extrapolation of L max /ro as a function of a/L max neglecting 
the error on the x-axis. We observe consistency between the two methods. 
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V. CONCLUSIONS AND DISCUSSIONS 

In this paper, we have calculated the step scaling function (SSF) at the weak and the strong couplings for both 
Iwasaki and LW actions with the one-loop 0(a) improved as well as the tree level 0(a) improved boundary terms. 
We have also calculated the low energy scale ratio for the Iwasaki action with both tree level and one- loop 0(a) 
improvements. The extrapolated values of the SSF at the weak and strong coupling for various gauge actions are 
consistent within la and 2.3a, respectively. The low energy scale ratio is also consistent between the Iwasaki and 
plaquette actions within la. In conclusion, we have confirmed the universality of both quantities. 

We have investigated lattice cut off effects in some detail. In the extrapolation procedure, the perturbative removal 
of lattice artifacts reduces the scaling violation of the SSF for the Iwasaki action with the tree level 0(a) improvement 
and the one- loop 0(a) improvement with the condition B. Indeed, at the strong coupling at the coarsest lattice 
L/a = 4, cut off effects are of order 1% and 3%, respectively, if one compares the extrapolated value obtained by the 
constrained fit. At the weak coupling, they are roughly 1% for both cases. We conclude that for the Iwasaki gauge 
action these combinations of improvements are the good choice for controlling lattice artifacts. This conclusion is 
also supported by the fact that an individual extrapolation to the continuum limit with the linear fitting form for the 
data set with the tree level 0(a) improvement or the one-loop 0(a) improvement with the condition B gives a result 
consistent with the extrapolated value estimated from the constrained fit within errors, at both weak and strong 
couplings. 

As mentioned in the introduction, this work is the second step toward Nf = 2 and 3 simulations. The present study 
shows that we should use the tree level 0(a) improved action or the one-loop 0(a) improved action with the condition 
B in future simulations with dynamical quarks. 
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FIG. 1. The assignments of the weight factor for loops near the boundary t = eq. (II. 5, II. 6). 
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FIG. 2. Results of SSF at the weak coupling u = 0.9944 with the Iwasaki action (upper) and LW action (lower). The 
simulation raw data (left-hand side) and the data with "perturbative removal of the lattice artifacts" (right-hand side) are 
shown. The points at a/L = in (b) and (d) represent the extrapolated values obtained by the constrained fit, that the dotted 
lines indicate the fitting function, for the Iwasaki action and by the linear fit, whose fitting function is shown as dotted line, 
for the data of 1-loop O(a) improved action for the LW action respectively. 
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FIG. 3. Results of SSF at the strong coupling u = 2.4484 with the various gauge actions. Concerning the symbols and the 
lines, the same explanation as the weak coupling case is followed. 
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FIG. 4. Interpolation of ro/a with phenomenological representation (eq. (IV. 7)) in the range 2.456 < f3 < 3.53. 
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on x-axis are invisible in this scale. The constrained fit is represented by dotted lines. In the fit the point at a/L max = 1/4 of 
the tree level O(a) improvement are not included. 
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action 


cf ^ with condition A 




cf ^ with condition B 


A 1 


Iwasaki 


0.1518 






0.04161 




0.3036 


LW 


-0.00297 










-0.00594 




TABLE I. The values of 


^ and Ai 


for the improved gauge actions with each condition. 




action 


degree of O(a) impr. 


L/a 


P 




9sp(2i) 


E(2,m,o/L) 




tree 


4 
6 
8 
12 


6.5447 
6.8485 
7.0733 
7.3765 


0.9944(5) 
0.9944(8) 
0.9944(12) 
0.9944(18) 


1.0953(14) 
1.0915(19) 
1.0973(27) 
1.0989(53) 


1.0953(16) 
1.0915(21) 
1 0973(311 
1.0989(57) 


Iwasaki 


one-loop 
condition A 


4 
6 
8 
12 


6.1467 

p. troop) 

6.8799 
7.2547 


0.9944(5) 

0.9944(13) 
0.9944(14) 


1.1395(16) 
1. IzoUI zl ) 
1.1192(29) 
1.1132(40) 


1.1395(17) 
1.1230(24) 
1 1192(331 
1.1132(43) 




one-loop 
condition B 


4 
6 
8 
12 


6.2258 
6.6358 
6.9010 
7.2722 


0.9944(5) 
0.9944(9) 
0.9944(9) 
0.9944(15) 


1.1284(15) 
1.1147(22) 
1.1133(30) 
1.1134(56) 


1.1284(16) 
1.1147(24) 
1.1133(31) 
1.1134(58) 


LW 


tree 


4 
6 
8 
12 


8.2189 
8.5889 
8.8479 
9.2017 


0.9944(4) 
0.9944(7) 
0.9944(10) 
0.9944(16) 


1.1177(12) 
1.1128(17) 
1.1136(27) 
1.1122(32) 


1.1177(13) 
1.1128(19) 
1.1136(30) 
1.1122(37) 


one-loop 
condition A 


4 
6 
8 
12 


8.2199 
8.5957 
8.8406 
9.2060 


0.9944(4) 
0.9944(7) 
0.9944(10) 
0.9944(16) 


1.1158(12) 
1.1115(17) 
1.1153(27) 
1.1085(47) 


1.1158(13) 
1.1115(19) 
1.1153(30) 
1.1085(51) 


TABLE II. Results of SSF at the weak coupling u = 0.9944. 
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action 


degree of O(a) impr. 


L/a 


13 


g 2 s AL) 


3sf(2L) 


S(2,m,o/L) 




tree 


4 
6 

8 

12 


3.2663 
3.5754 
3.7872 
4.0996 


2.4484(28) 
2.4484(52) 
2 4484(40 ~) 

. 1 IV.,' J I IV.' / 

2.4484(65) 


3.332(11) 
3.352(16) 
3 361(24^ 
3.396(43) 


3.332(12) 
3.352(18) 
3 361 f 251 
3.396(44) 


Iwasaki 


one-loop 
condition A 


4 
6 

8 

12 


2.9628 
3.3803 
3.6544 
4.0091 


2.4484(29) 
2.4484(60) 
2.4484(42) 
2.4484(68) 


4.008(21) 
3.692(21) 
3 625('27~) 
3.512(45) 


4.008(21) 
3.692(23) 
3 625(281 
3.512(46) 




one-loop 
condition B 


4 
6 
8 
12 


3.0624 
3.4395 
3.6908 
4.0283 


2.4484(29) 
2.4484(56) 
2.4484(47) 
2.4484(64) 


3.712(15) 
3.573(20) 
3.519(28) 
3.518(38) 


3.712(16) 
3.573(21) 
3.519(29) 
3.518(39) 


LW 


tree 


4 
6 
8 
12 


4.8992 
5.2786 
5.5325 
5.8878 


2.4484(24) 
2.4484(49) 
2.4484(42) 
2.4484(62) 


3.619(13) 
3.540(18) 
3.486(23) 
3.468(32) 


3.619(14) 
3.540(20) 
3.486(24) 
3.468(33) 


one-loop 
condition A 


4 
6 
8 
12 


4.9055 
5.2784 
5.5332 
5.8867 


2.4484(25) 
2.4484(51) 
2.4484(40) 
2.4484(58) 


3.601(13) 
3.530(18) 
3.499(26) 
3.463(25) 


3.601(13) 
3.530(19) 
3.499(26) 
3.463(27) 




TABLE III. 


Results of SSF at the strong coupling u = 2.4484. 







degree of O(a) impr. 




P 


5SF(^max) 


tree 


4 
6 
8 


2.7000 
3.0057 
3.2154 


3.480(6) 

3.480(11) 

3.480(11) 




12 


3.5219 


3.480(13) 




4 


2.4594 


3.480(7) 


one-loop 


6 


2.8556 


3.480(14) 


condition A 


8 


3.1047 


3.480(11) 




12 


3.4496 


3.480(16) 




4 


2.5382 


3.480(6) 


one-loop 


6 


2.8921 


3.480(11) 


condition B 


8 


3.1376 


3.480(12) 




12 


3.4734 


3.480(14) 



TABLE IV. Tuning of f3 at u = 3.480 for Iwasaki action. 



I (L/a) 4 ro/a Nor iVconf 



3.00 


32 4 


8.88(13) 


15 


400 


3.53 


56 4 


17.35(13) 


24 


88 



TABLE V. Simulation parameters and results performed in this work for ro/a with the Iwasaki action. A^or and A con f 
indicate the number of over-relaxation sweep and configuration respectively. 
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p 


ro/a 


reference 


2.456 


4.080(16) 


[29] 


2.461 


4.089(14) 


[28] 


2.487 


4.286(15) 


[29] 


9 KOQ 
Z.OZo 


4.0 (KJ{Z± ) 


[29J 


2.575 


4.887(16) 


[29] 


2.659 


5.556(30) 


[28] 


3.000 


8.88(13) 


in this work 


3.200 


11.53(15) 


[27] 


3.530 


17.35(13) 


in this work 


TABLE VI. 


(3 versus ro/a with the Iwasaki action. The value of ro/a are 


taken from the references quoted in the last column. 





Iwasaki action 




LW action 




L/a 


5 (1A) 








4 -0.02096 
6 -0.01922 
8 -0.01499 
12 -0.01064 


0.01700 
0.00608 
0.00399 
0.00201 


0.01577 
0.00592 
0.00395 
0.00200 


0.003278 
0.000911 
0.000527 
0.000296 


0.002536 
0.000417 
0.000156 
0.000049 



TABLE VII. 1-loop deviations for various gauge actions and with fc-level O(a) improvement. 



action 


a(2,u = 0.9944) 




a(2,u = 2.4484) 


xV^dof 


Iwasaki 


1.107(3) 


1.5/8 


3.485(34) 


1.9/6 


LW 


1.111(4) 


2.0/2 


3.410(30) 


0.1/2 


plaquette [17] 


1.110(11) 


3.3/2 


3.464(40) 


0.9/4 


TABLE VIII. The 


extrapolated values of the SSF with various gauge actions for both the weak and strong coupling. The 


values in the case of the plaquette action are quoted from the reference as 


indicated. 




degree of O(a) impr. 


a 


r /a 


^max / a 


imax/ ^0 




2.7000 


5.886(26) 


4.000(11) 


0.680(4) 


tree 


3.0057 


9.03(12) 


6.000(35) 


0.664(10) 




3.2154 


11.87(19) 


8.000(41) 


0.674(12) 




3.5219 


17.16(14) 


12.00(8) 


0.699(7) 




2.4594 


4.098(12) 


4.000(9) 


0.976(4) 


one-loop 


2.8556 


7.352(59) 


6.000(33) 


0.816(8) 


condition A 


3.1047 


10.29(16) 


8.000(31) 


0.777(12) 




3.4496 


15.78(16) 


12.00(7) 


0.761(9) 




2.5382 


4.625(12) 


4.000(8) 


0.865(3) 


one-loop 


2.8921 


7.735(73) 


6.000(28) 


0.776(8) 


condition B 


3.1376 


10.74(17) 


8.000(38) 


0.745(12) 




3.4734 


16.22(15) 


12.00(7) 


0.740(8) 



TABLE IX. The low energy scale at various values of /3 for Iwasaki action. 
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action 




xV^dof 


Iwasaki 


0.749(14) 


3.93/ 5 


plaquette [18] 


0.738(16) 





TABLE X. The extrapolated values of the low energy scale ratio with various gauge actions. The values in the case of the 
plaquette action are quoted from the reference as indicated. 
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